Fluctuation-Dissipation relations in Driven Granular Gases 



A. Puglisi, A. Baldassarri and V. Loreto 
La Sapienza" University in Rome, Physics Department, P.le A. Moro 5, 00185 Rome, Italy and 
INFM, Center for Statistical Mechanics and Complexity, Rome, Italy 
(Dated: February 1, 2008) 

We study the dynamics of a 2d driven inelastic gas, by means of Direct Simulation Monte Carlo 
(DSMC) techniques, i.e. under the assumption of Molecular Chaos. Under the effect of a uni- 
form stochastic driving in the form of a white noise plus a friction term, the gas is kept in a 
non-equilibrium Steady State characterized by fractal density correlations and non-Gaussian distri- 
butions of velocities; the mean squared velocity, that is the so-called granular temperature, is lower 
than the bath temperature. We observe that a modified form of the Kubo relation, which relates the 
autocorrelation and the linear response for the dynamics of a system at equilibrium, still holds for 
the off-equilibrium, though stationary, dynamics of the systems under investigation. Interestingly, 
the only needed modification to the equilibrium Kubo relation is the replacement of the equilibrium 
temperature with an effective temperature, which results equal to the global granular temperature. 
We present two independent numerical experiment, i.e. two different observables are studied: (a) 
the staggered density current, whose response to an impulsive shear is proportional to its autocorre- 
lation in the unperturbed system and (b) the response of a tracer to a small constant force, switched 
on at time t m , which is proportional to the mean-square displacement in the unperturbed system. 
Both measures confirm the validity of Kubo's formula, provided that the granular temperature is 
used as the proportionality factor between response and autocorrelation, at least for not too large 
inelasticities. 

PACS numbers: 05.40.-a,45.70.-n,05.70.-a 



I. INTRODUCTION 

In the past few years granular materials jl| have become a fast growing field of research. In this framework, a 
strong interest has arisen in the subject of the so-called granular gases, systems of grains at very low density so that 
the collisions can be considered always binary and instantaneous Q . The collisions among the grains are dissipative, 
and the amount of energy lost in each collision is parametrized by the so-called restitution coefficient r (see below for 
a precise definition). Since the energy of the system is therefore not conserved an external forcing is usually applied 
to obtain a dynamic stationary state. For this kind of systems, the analogy with perfect gases allows to introduce a 
granular temperature || defined as the mean square fluctuation of the velocity, that is Tq = (|v — (v)| 2 ) even though 
the velocity distribution could not be Gaussian. 

In the case of zero external forcing [5J (i.e. the so-called free cooling), granular temperature decreases in time and 
goes asymptotically to zero when all the particles finally stop. When an energy input feeds the system, instead, the 
granular temperature may reach a stationary value and several modeling for the stochastic driving have been proposed 
employing a constant || ^, [j], |, | or a random restitution coefficient. 

In numerical simulations or in experiments, Tq is often measured taking the average < • > on the whole system. 
However strong fluctuations in the local granular temperature are usually observed Q. Even considering only the 
global granular temperature, i.e. for small spatial inhomogeneities, it is not clear to which extent Tq can be considered 
the "temperature" of the system. 

There are of course several possible paths to face this problem. One interesting point of view is that of investigating 
the response properties of an external thermometer coupled to a granular gas pd| ]. This corresponds to study the 
fluctuation-dissipation properties of the system fl^| . For a system slightly perturbed from its stationary equilibrium 
state linear response theory allows to relate the response to the correlation functions through the fluctuation-dissipation 
relations Jlq] . 

In the simplest case, given a perturbing field a, a fluctuation-dissipation relation relates the response of an observable 
B at the time t, after an impulsive perturbation at time 0, to the correlation of the observable B and the field A, 
conjugated to a, measured in the unperturbed system, 

d(B(t)) 19 

where T is the equilibrium temperature of the system. 

Recently Green-Kubo expressions for a homogeneous cooling granular gas have been obtained jL4|, |l6| and 
numerical simulations have confirmed their validity J17[: cooling granulars lack an equilibrium state, therefore the 
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Homogeneous Cooling State (which is characterized by scaling properties) is used as reference state to be perturbed, 
but in this case the Green-Kubo relations must be changed in order to keep into account new terms arising from the 
time dependence of the reference state and the non-conservative character of collisions (see Jl5|). However, in the 
case of steady state granular gases, a rigorous derivation of Green-Kubo relations has not yet been performed to our 
knowledge. 

In this paper we perform numerical investigations in order to check the validity of standard Kubo formulas | fl3[ |l§f| 
to steady state inelastic gases. We perform in particular two different sets of numerical experiments on heated granular 
gases. We choose two different conjugated pairs of variables constituted by the autocorrelation of a given variable and 
the corresponding response to a perturbation applied to the gas. Kubo's formula, i.e. the proportionality between 
response and the autocorrelation in the unperturbed system, is verified to hold using the granular temperature as the 
correct proportionality factor, at least for not too strong inelasticities (r > 0.5). For stronger inelasticities (smaller 
values of r) Kubo's formula can be verified with a different effective temperature. It should be noted how the adopted 
simulation scheme (i.e. Direct Simulation Monte Carlo) represents the numerical implementation of the Boltzmann 
equation (given for example in [Q) which assumes Molecular Chaos hypothesis and therefore neglects short range 
correlations. For this reason, for very strong inelasticities the Boltzmann equation (and the DSMC scheme) cannot 
be considered realistic and one should use Molecular Dynamics simulations. 

These results differ from the analogous results obtained for dense granular assemblies, where "slow" degrees of 
freedom thermalize at an effective temperature which is far higher than the external imposed temperature p0[ . 
However the gas-like state of granular matter has nothing to do with these systems, as their stationary state is 
governed by a rapid decay of the fluctuations and the granular temperature turns out to be the right choice in the 
description of linear response to slight perturbations. 

The outline of the paper is as follows. In section II we define the model as well as its simulation scheme and give a 
brief review of known results about the peculiar features of its non-equilibrium stationary state. Section III reviews 
Kubo's formula and its extension to non-Hamiltonian perturbations. In section IV we describe numerical experiments 
using a sinusoidal shear and measuring as response the density current. In section V we perform a diffusion vs. 
mobility experiment upon a tracer (i.e. perturbing a single particle). Finally section VI is devoted to the conclusions. 



II. THE MODEL 



We simulate a gas of N identical particles of unitary mass in a two-dimensional box of side L — \/~N with periodic 
boundary conditions. The particle collisions are inelastic: the total momentum is conserved, while the component 
of the relative velocity parallel to the direction joining the center of the particles is reduced to a fraction r (with 
< r < 1) of its initial value, lowering in this way the kinetic energy of the pair: 



1 + r 

v'i=vi — ((vi - v 2 ) -n)n (2a) 

V2=v 2 + ^i^(( Vl -v 2 ).n)n (2b) 

where and v 2 are the velocities of the colliding particles after the collision. 

In the interval between two subsequent collisions, the motion of each particle i is governed by the following Langevin 
equation |t]]: 



— Vi(t) = + J tjM 3a 

at Tb V TB 

^x,(t)=Vi(t) (3b) 

where the function rj^t) is a stochastic process with average (fj^t)) — and correlations (r)f(t)r^ (t')) = S(t—t')Sij5 a p 
(a and /3 being component indexes), i.e. a standard white noise. This means that each particle feels a hot fluid with 
a temperature Tb and a viscosity characterized by a time tb- 

The question about the most proper way of modeling a stochastic driving is still open. Many authors for instance 
use Eq. (||) without viscosity M. This is equivalent to the limit tb — ► oo and Tb — ► oo with keeping D = Tb/tb 
constant. In this limit long and short-range correlations in the velocity and density fields have been observed B. We 
have measured correlations in the velocity field in the model with viscosity, concluding that they are highly reduced 
by the viscous term that breaks Galilean invariance (the frame v = is preferred) p3 . Our choice of the stochastic 
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driving with viscosity has the following advantages: a) it guarantees that in the elastic limit (r = 1) the system, 
after a transient time of the order of Tg, still reaches a stationary state, characterized by a uniform density and a 
Gaussian distribution of velocities whith temperature Tb; b) it is a heat bath with a well defined finite temperature 
Tb that can be compared with the granular temperature Tq and the effective temperature T e f f measured by means 
of Fluctuation-Dissipation relations (see ahead) . 




FIG. 1: Rescaled (in order to have unitary variance) distribution of the horizontal component of the velocity of the gas P(v x ) 
versus v x , with N — 500, L 2 — N, heat bath with Tb = 1, Tb = 10, r c = 1, and different restitution coefficients. The Gaussian 
is plot as a reference for the eye. 



In spite of the drastic reduction of velocity correlations, this model exhibits several interesting features which are 
parametrized by the restitution coefficient r and the ratio, tc/tb, between the mean free time (the average interval 
of time between two subsequent collisions of the same particle) and the viscosity time. When tb > tc and r < 1 the 
system is in a non-equilibrium stationary state with a granular temperature Tq < Tb- This state is characterized R] 
by fractal density clustering and non-Gaussian distributions of velocities (see fig.([lj)) with a dependence upon r and 
tc/tb- This situation persists when Direct Simulation Monte Carlo (DSMC) |22| is used to perform more rapid and 
larger simulations of the system. The DSMC simulation scheme consists of a discrete time integration of the motion 
of the particles. At each time step of length At the following operations are performed: 

1. Free streaming: Eq. is integrated for a time step At disregarding possible interactions among the particles. 

2. Collisions: every particle has a probability p — At/r c of undergoing a collision. Its collision mate is chosen 
among the particles in a circle of fixed radius tb with a probability proportional to its relative velocity. The 
unitary vector n, which should be parallel to the line joining the centers of particles, is chosen randomly instead 
and the collision is performed using rule (0). 



III. SHORT REVIEW ABOUT KUBO'S FORMULA 



In this section we review the Kubo's formula, which must hold in the case of equilibrium dynamics. First we show 
the case of an Hamiltonian dynamics, and then for non-Hamiltonian equilibrium dynamics of a system subjected to 
an (impulsive) shear force. The first case is presented just for completeness and to set the notations, the formulas 
presented in the second case will be numerically checked for an elastic system and extended to the inelastic case in 
the next sections. 



A. Kubo's formula for Hamiltonian systems 



For an Hamiltonian system, the temporal evolution of the phase space distribution f(p, q, t) is ruled by the Liouville 
equation: 



-gif(p,q,t) = 



£ext(t)]f(p, q,t) 



(4) 
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where C and C ex t are the Liouville operators relative to the unperturbed Hamiltonian and to its perturbation respec- 
tively. They are defined by means of classical Poisson bracket: 

The perturbation is assumed to be given by coupling a force a(t) with an observable of the system A, i.e.: 

n ex t = -A(p,q)a(t) (6) 

The phase space distribution is assumed to be canonical (i.e. at equilibrium) in the infinite past: f(p, q, — oo) = 
f e q(p, ?)• An approximate solution of (Q) to the first order in the perturbation i£ ex t is given by: 

f(p,q,t)=f eq (p,q)+ I dt'e^[i{t-t')L]{H ext (t'),f eq (p,q)) + ... (7) 



With this approximation, the deviation, due to the perturbation (g), of the expectation of a physical quantity B 
can be written as the following convolution: 

SB(t) = (B) t - (-B)-co = / dt'$ BA {t - t')a{t) (8) 



where (■ ■ - )t denotes averages taken over the ensemble given by f(p,q,t); $BA(t — t'), called response function, 
represents the response SB(t) to the pulsed force a(t) = S(t), and reads: 

$BA(t) = J J dpdqf eq (p,q)(AA(p,q),AB(p t ,q f )) = ((AA,AB(t))) (9) 

where A A = A — (A) (and identically for AB), while (pt, qt) is the image of the initial phase point (p, q) determined 
by the total Hamiltonian (including the perturbation). 

Kubo has shown that the response function can be written in a simpler form [l^] : 

f) A A f)/\ H 

*ba® = ^(-g r (0)AB(t)) = -/?<AA(0)— (t)) (10) 

where (3 is the inverse of the temperature. 

B. Fluctuation-Dissipation for non-Hamiltonian equilibrium systems (the case of shear force) 

Let consider a gas of particles and define a non-conservative perturbation (force) acting on particle i placed at (t) 
at time t as 

F(r i; t) = 7i£(rj, t) with the properties Vx^O V ■ £ = 0. (11) 

For small enough perturbation and for any variable ("observable") 0(r) such that (O)-oo = 0, equation (^|) with the 
Kubo formula @ reads ||, |§: 

(0(r)>« =0 [dr> [ dt'{0{v, t) ]T nitfW - ■ £(r', t') (12a) 

(6(k)) t =/3 f dt , <6(k,t)^ 7i r i (t')exp{ikT i (t')})-oo^(k,t') (12b) 

J —OO 
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where the Fourier transform G — > G is defined as 



G(k) = ^ J dve-^G{v) 
A force satisfying the properties ((O)) is for instance given by: 



(13) 





S exp(ifc x a;)(5(t) 



(14) 



whose spatial Fourier transform reads: 



€s(M) 



o _ 

^(k-k)tf(t) 



(15) 



where k = (fc x ,0) (having chosen k x compatible with the periodic boundary conditions, i.e. k x = 2nkTt/L x with rik 
integer and L x the linear horizontal dimension of the box where the particles move). Note that S must have the 
dimensions of a momentum, i.e. of a velocity (taking unitary masses). With this choice, eqs. (12t) becomes: 



(0(k)) ( = ^<0(k, t) 7iW(0) exp{zk • r i (0)})_ oo «(k - k) 

i 

If we now define the staggered y-current as: 

4 t (r,t) = ^ 7i y i (^(r-r i (t)) 

i 

J*(k,f) = 1^7 i y i (i)exp(-tkT i (t)), 

i 

then, using this current as observable O the relation (|l6| ) is written as: 



(16) 

(17a) 
(17b) 



«(M))t 



r 2 



C^liVi{t) exp(-ik ■ r i (t))7 J ?) J -(0) exp{i/c :E x- ) -(0)})_ oo (5(k - k) 



(18) 



= /3s( r y l (k, t) (-k, o)) _«,*(]£ - k) 

A real linear combination of forces of the kind in eq. ( |l4| ) is 

C(r,t) = ~(fc(r,f)+*_ £ (r,t)) = (g^^ 
With this choice of the perturbation, the relation ( |l8| ) becomes: 

(jf (k,t)) ( = ^<i y rt (k,t)j«(-k,o))_ 00 («(k-k)+*(k+k)). 



(19) 



(20) 



This is a fluctuation- dissipation relation which expresses the fact that the response of the k-component of the transverse 
current to the perturbing field in eq. ( |l9| ) is proportional to the auto-correlation of the same transverse current 
measured in the system without perturbation. 

The real part of the response calculated at k (per unit of perturbing field) is directly computable and reads: 



Re 



V 



= — (^Ti&C*) cos(k x Xi(t))) t . 



From eq.(20) one obtains the relation: 



^C^ltMt) cos(k xXl (t))) t = ^(^2i t ijm(t)yj(0) cos{k x [x t (t) - a!j(0)]}>. 



(21) 



(22) 
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IV. FLUCTUATION-DISSIPATION MEASURE I: SHEAR VS. CURRENT 

The first set of measures we have performed has been aimed to verify relation ( |22| ) for a system as described in the 
previous paragraph. In order to do this we have used the following recipe proposed in p3[ : 

1. Initialize system U with random positions {^(0)} and random velocities {^(O)}. 

2. Let it evolve with the unperturbed dynamics until time t w which must be chosen larger than the largest 
characteristic time of the system (e.g. t c or r^). The unperturbed dynamics consists of the time-discretized 
(At) integration of the Langevin equation (||): 



Vi(t + At) = Vi (t) - — v«(t) + «p^R(t) (23) 
n V n 

plus inelastic collisions with parameter r (restitution coefficient). The collisional step is separated from the 
Langevin step and is implemented by means of a local Monte Carlo, i.e. random choice of pairs to collide inside 
a region of diameter re; the collision probability is proportional to a fixed collision frequency l/r c and to the 
relative velocity of the particles; r c is chosen to be compatible with an homogeneous gas-like dynamics, i.e. 

T c ~ r B f '\f{v 2 )-oo- 

3. At time t w a copy of system U is created (and named P) and the vectors {yi(t w )} and {xi(t w )} memorized in 
order to be used in the computation of the auto-correlation. 

4. The system U is let evolve with the unperturbed dynamics. The system P is made evolve with the additional 



forcing described in eqs. (11) and (113) for only the time step [t w ,t w + At], i.e. the equation for the update of 



velocities in this particular step being: 



vf (t w + At) = vf (t w ) - ^vf (t w ) + f^m w ) + H [ Ecos[ k° xX ? (t w )}) ^ 

with k x = 2irnk/L x . Note again that the perturbation intensity S has exactly the dimensions of a velocity. 

5. The dynamics of systems U and P are thereafter followed in the unperturbed style, i.e. using eq. (|23|). The 
functions to be measured are: 



R{t) = -j^=Yl ^ ( tw + T ) cos ^ x i + T )) ( 25a ) 

i 

C ( r ) = jf 12™?^ + T)yJ(t™) cos{k x [xY(t w + t)- xf (t w )}}, (25b) 

where T = t-t w . It is expected that i?(0) = 1/2 and C(0) = (v 2 ), while R(oo) = C(oo) —> 0. 

6. The above steps are repeated for many different realizations (or even in the same realization, provided its length 
is much longer than the typical correlation time) and the averages over those realizations of R(t — t w ) and 
C(t — t w ) are computed (see figure ||). 

To the aim of checking the whole numerical machinery, we first consider the r = 1 elastic case. In all the cases 
investigated we have checked the linearity of the response by changing the perturbation amplitude in the range 
H € [0.005,0.05]. Within the specific framework chosen for the observables, the Kubo formula to be verified is given 
by: 



(R(t-t w )) = ^(C(t-t w )) (26) 



We have performed the following experiments: 
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FIG. 2: Left: time dependent response to the impulsive shear perturbation defined in [ll] and |l^: R(t — t w ) vs. t — t w for three 
simulations with elastic systems, one without thermal bath and two with thermal bath, and with different choices of the wave 
number of the perturbation. Right: time correlation function C(t — t w ) vs. t — t w for the same systems. In the simulations 
r = 1, N = 500, r c = 1, and for the two cases with the heat bath Tb — 1 and Tb = 10. The applied force has H = 0.01, Tw = 100. 
The averages have been obtained over 10000 realizations. 




C(H„ ) 



FIG. 3: Parametric plot of R(t — t w ) vs. C(t — t w ) for the numerical experiment of type I (impulsive shear perturbation) 
with r = 1, with or without heating bath, and for different choices of the wave number of the perturbation. The initial 
temperature for the case without heat bath is chosen to be 1 while T;, = 1 ad rj = 10 for the two cases with the heat bath. 
N — 500, t c = 1, H = 0.01, rik = 8, with average over 10000 realizations, using t w = 100. 

• gas with elastic interactions and absence of thermal bath (see figure ||) . 

• gas with elastic interactions with the thermal bath (see figure ||) . 

From figure |^ it can already be appreciated that response and autocorrelations in the elastic gas decay on a time of 
the order of r c (however this decay is not exponential, as can be observed in the inset of the figure) . The parametric 
plot of the two curves in figure [| shows the perfect agreement with equation ( p6| ) using (3 = 1/X&. In this case, of 
course, Tb = Tq. 

We have repeated the same measures on the gas with restitution coefficient r < 1, i.e. 

• gas with inelastic interactions with the bath (in this case the bath is essential, to avoid cooling) , see figure [I| 

obtaining again a very good agreement with equation (|2^) using (3 — \/Tq. This is the main finding in this set 
of numerical experiments: even if the gas is out of equilibrium, being driven by a thermal bath at temperature Tb, 
its unperturbed autocorrelation is still proportional to the linear response, and its effective temperature measured by 
means of Fluctuation-Dissipation theory is exactly the granular temperature Tq. This is true for every inelasticity 
r > 0.5. 

At lower inelasticities we have still obtained an agreement with equation (^) but with a T e f f = 1//3 ^ Tq : we do 
not think that these measures can be interesting, as the DSMC algorithm is not reliable at such values of r. A real 
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system of inelastic particles with r < 0.5 should exhibit, even at low densities, a high degree of clusterization which 
cannot be observed within the DSMC framework. 

It also must be noted that relation ( p6| ) is verified for many values of the wave number rife, i.e. the system does not 
show a scale dependent effective temperature. 
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FIG. 4: Parametric plot of R(t — t m ) vs. C(t — t w ) for the numerical experiment of type I (impulsive shear perturbation) with 
r < 1, with heating bath, and for different choices of the wave number n k of the perturbation. Tj, = 1 and Tt = 10, N = 500, 
t c = 1, S = 0.01, n k — 8, with averages over 10000 realizations, using t w = 100. 



V. FLUCTUATION-DISSIPATION MEASURE II: DIFFUSION 



Another independent confirmation of the validity of the modified Linear Response theory for granular gas comes 
from the study of the diffusion properties, i.e. of the large time behavior of the mean squared displacement B(t, t w ) = 
(|r(<) — r(t w )\ 2 ) ~ 2D(t — t w ). In this case the Einstein relation is expected to hold D =< |v| 2 > T corr , where 
Tcorr = ft J dr{v{t w + r)v(t w )}: this relation however, often addressed as a sort of FD relation, is always verified 
and just represents a check of the correctness of the simulation. Instead some surprise could arise from mobility 
measurements: a small static drag force (switched on at time t w ) of intensity 3 (in the direction of the unitary 
vector x of the x axis) is applied to a given particle (tracer, e.g. particle with index and position r ). The tracer 
reaches, as a result of the viscous force generated by the gas surrounding it, a limit constant velocity such that 
x(t,t w ) = (|(ro(i) — vo(t w )) ■ x\) ~ S/ii, where n is the mobility which is expected to be related to the diffusion 
coefficient through the Einstein relation /i = D/(v 2 ) = 2D/T (if the force is applied on the direction x in the 
two-dimensional system). 

In our experiments we have checked the linearity of the relation between (|^o(*) — %o{t w )\) = (|( r o(£) — ro(iw)) ' ^1) 
and (|r(t) — r(t w )\ 2 } (see fig. ||). In particular if Kubo's formula was valid one should have: 



(\x (t) -x (t w )\) = p (\rj(t) - r^l 2 ) (2?) 

with = 1/T if the system is in thermodynamic equilibrium at temperature T. 

In all the simulations we have checked the linearity of the response by changing the perturbation amplitude in the 
range £ €[0.005,0.05]. 

Figure | shows the mean squared displacement (in the unperturbed system) and the x displacement of the tracer 
(when it is accelerated) divided by the intensity of the perturbing force versus time: it can be appreciated how both 
these quantities grow linearly with time, defining the diffusion coefficient and the mobility. 

Figures (^,0) report the parametric plots of the response to the force versus the mean squared displacement in the 
unperturbed system, showing how relation ( p7| ) is very well satisfied at different inelasticities. What can be observed 
in this kind of measures is a departure from the relation (3 = 1/Tq already at inelasticity r = 0.5 (i.e. where the 
check of fluctuation-dissipation theorem with the first type of measures still was positive): it appears that (3 = 1/T e ff 
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with T;, > T e tf > Tq- This breakdown of fluctuation-dissipation relation should be related to cluster formation (i.e. 
spatial lack of homogeneity), which is present even in the case of DSMC solutions of Boltzmann equation (see for 
example 0). This violation is more pronounced than in the previous experiment, therefore appearing earlier, because 
this is a "local" measurement of FD, i.e. only one particle is perturbed: this means that the fluctuations of the local 
granular temperature strongly influence the trajectory of the tracer and its statistical properties. 




20 40 60 80 100 20 40 60 80 100 
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FIG. 5: Simulations of systems coupled to a thermal bath with elastic or inelastic collisions. TV = 500, r c = 1, Tj = 0.1 and 
Tt = 10, S = 0.01, t w — 100. The results are obtained by averaging over 10000 realizations. Left: mean squared displacement 
B(t,t w ) vs. t — tw Right: Integrated response x(*>£™) to a constant force applied to the particle numbered as vs. t — t w . 




FIG. 6: Parametric plot of x(t> tw) vs. B(t, t w ) for the numerical experiment of type (constant force applied on one particle) 
with r = 1, r = 0.8 and r — 0.7, with heating bath, and for different choices of the intensity H of the perturbation, using 
Tb = 0.1 ad Tb = 10, N = 500, r c = 1, t w = 100. The results are obtained by averaging over 10000 realizations. 



VI. CONCLUSIONS 



In conclusion in this paper we have investigated the validity of Kubo's relations in Driven Granular Gases in 
two dimensions. We have compared in particular two sets of measures. On the one hand we have measured the 
proportionality factor between the response of the staggered density current to an impulsive shear forcing and its 
autocorrelation function in the unperturbed system. On the other hand we have monitored the velocity of a tracer 
by checking the proportionality between its response to a small constant force, switched on at time t w , and its mean 
squared displacement in the unperturbed stationary state. 

In both cases a proportionality is observed, in analogy with the Linear Response theory for equilibrium dynamics. 
Furthermore the proportionality factor in the Kubo formulas is equal to the inverse of th granular temperature (at 
least for the restitution coefficients larger than r > 0.5), which hence plays the role of the equilibrium temperature in 
the elastic case. It is important to remark how these results are recovered by two completely independent measurement 
schemes. 
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FIG. 7: Parametric plot of x(t> t w ) vs. B(t, t w ) for the numerical experiment of type (constant force applied on one particle) 
with r — 0.5, with heating bath, and for different choices of the intensity E of the perturbation, using Tb = 0.1 ad Tb = 10, 
N = 500, t c = 1, with average on 10000 realizations, using t w = 100. This case r = 0.5 exhibits a slope compatible with 
T eff = 0.051 which is above Tc- 



Several remarks are in order. First of all, though a granular gas is a non-trivial out-of-equilibrium system, from 
the point of view of its thermodynamics it exhibits properties which seem much simpler than the corresponding 
properties observed in a compacting granular medium. In this case in fact apparently no slow modes are present or at 
least their presence does not give raise to the existence of several effective temperatures depending on the time-scales 
investigated JO ], On the other hand the proportionality factor between response and autocorrelation we have found 
cannot be considered a temperature from the point of view of equilibrium thermodynamics since it does not rely on 
any known statistical ensemble. 

Another point to stress concerns the validity of the zeroth principle of thermodynamics. The question that imme- 
diately arises could be summarized as follows: if the granular temperature represents the correct temperature from 
the point of view of the Fluctuation-Dissipation Theorem, should we expect it rules the thermalization properties 
of two different granular gases put in contact? As already mentioned the answer to this question is far from being 



trivial as witnessed by all recent results obtained on mixtures [|25|, |26| |27j, |28|, [29|, |30|, |31|, |32j, |33|, |3J, |35J where a lack of 
equipartition is observed. However these results (validity of fluctuation-dissipation relations and lack of equipartition) 
can coexist simply because heated granular gases have not only a thermal source but also a thermal sink (dissipative 
collisions) and therefore any zero principle should be stated in terms of a balance equation among energy fluxes instead 
of a strict equivalence between temperature of systems in contact. 

We again remark, also, that there are other common ways of driving a granular gas in a stationary state, e.g. 
stochastic driving without viscosity ||, stochastic restitution coefficient models [[to), multiplicative noise models || 
and so on. In some of these models a more pronounced departure from homogeneity (for example correlations in the 
velocity field B) and therefore a breakdown of fluctuation-dissipation relations should be investigated. 
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